On the stability of the atmosphere-vegetation system in the Sahara/Sahel region

A conceptual model has been developed for the analysis of atmosphere-vegetation interaction in subtropical deserts. The model can exhibit multiple stable states-in the system: a "desert" equilibrium with low precipitation and absent vegetation and a "green" equilibrium with moderate precipitation and permanent vegetation cover. The conceptual model is applied to interpret the results of two climate-vegetation models: a comprehensive coupled atmosphere-biome model and a simple hox model. In both applications, two stable states exist for the western Sahara/Sahel region for the present-day climate, and the only green equilibrium is found for the mid-Holocene climate. The latter agrees well with paleoreconstructions of Sahara/Sahel climate and vegetation. It is shown that for present-day climate the green equilibrium is less probable than the desert equilibrium, and this explains the existence of the Sahara desert as it is today. The difference in albedo between the desert and vegetation cover appears to be the main parameter that controls an existence of multiple stable states. The Charney's mechanism of self-stabilization of subtropical deserts is generalized by accounting for atmospheric hydrology, the heat and moisture exchange at the side boundaries, and taking into account the dynamic properties of the surface. The generalized mechanism explains the self-stabilization of both desert and vegetation in the western Sahara/Sahel region, The role of surface roughness in climate-vegetation interaction is shown to be of secondary importance in comparison with albedo. Furthermore, for the high albedo, precipitation increases with increasing roughness while, for the low albedo, the opposite is found

A simple conceptual model of abrupt glacial climate events

Here we use a very simple conceptual model in an attempt to reduce essential parts of the complex nonlinearity of abrupt glacial climate changes (the so-called Dansgaard-Oeschger events) to a few simple principles, namely (i) a threshold process, (ii) an overshooting in the stability of the system and (iii) a millennial-scale relaxation. By comparison with a so-called Earth system model of intermediate complexity (CLIMBER-2), in which the events represent oscillations between two climate states corresponding to two fundamentally different modes of deep-water formation in the North Atlantic, we demonstrate that the conceptual model captures fundamental aspects of the nonlinearity of the events in that model. We use the conceptual model in order to reproduce and reanalyse nonlinear resonance mechanisms that were already suggested in order to explain the characteristic time scale of Dansgaard-Oeschger events. In doing so we identify a new form of stochastic resonance (i.e. an overshooting stochastic resonance) and provide the first explicitly reported manifestation of ghost resonance in a geosystem, i.e. of a mechanism which could be relevant for other systems with thresholds and with multiple states of operation. Our work enables us to explicitly simulate realistic probability measures of Dansgaard-Oeschger events (e.g. waiting time distributions, which are a prerequisite for statistical analyses on the regularity of the events by means of Monte-Carlo simulations). We thus think that our study is an important advance in order to develop more adequate methods to test the statistical significance and the origin of the proposed glacial 1470-year climate cycle

The role of orbital forcing, carbon dioxide and regolith in 100 kyr glacial cycles

The origin of the 100 kyr cyclicity, which dominates ice volume variations and other climate records over the past million years, remains debatable. Here, using a comprehensive Earth system model of intermediate complexity, we demonstrate that both strong 100 kyr periodicity in the ice volume variations and the timing of glacial terminations during past 800 kyr can be successfully simulated as direct, strongly nonlinear responses of the climate-cryosphere system to orbital forcing alone, if the atmospheric CO<sub>2</sub> concentration stays below its typical interglacial value. The existence of long glacial cycles is primarily attributed to the North American ice sheet and requires the presence of a large continental area with exposed rocks. We show that the sharp, 100 kyr peak in the power spectrum of ice volume results from the long glacial cycles being synchronized with the Earth's orbital eccentricity. Although 100 kyr cyclicity can be simulated with a constant CO<sub>2</sub> concentration, temporal variability in the CO<sub>2</sub> concentration plays an important role in the amplification of the 100 kyr cycles

Toward generalized Milankovitch theory (GMT)

In recent decades, numerous paleoclimate records and results of model simulations have provided strong support for the astronomical theory of Quaternary glacial cycles formulated in its modern form by Milutin Milankovitch. At the same time, new findings have revealed that the classical Milankovitch theory is unable to explain a number of important facts, such as the change in the dominant periodicity of glacial cycles from 41 to 100 kyr about 1 million years ago. This transition was also accompanied by an increase in the amplitude and asymmetry of the glacial cycles. Here, based on the results of a hierarchy of models and data analysis, a framework of the extended (generalized) version of the Milankovitch theory is presented. To illustrate the main elements of this theory, a simple conceptual model of glacial cycles was developed using the results of an Earth system model, CLIMBER-2. This conceptual model explicitly assumes the multistability of the climate–cryosphere system and the instability of the “supercritical” ice sheets. Using this model, it is shown that Quaternary glacial cycles can be successfully reproduced as the strongly nonlinear response of the Earth system to the orbital forcing, where 100 kyr cyclicity originates from the phase locking of the precession and obliquity-forced glacial cycles to the corresponding eccentricity cycle. The eccentricity influences glacial cycles solely through its amplitude modulation of the precession component of orbital forcing, while the long timescale of the late Quaternary glacial cycles is determined by the time required for ice sheets to reach their critical size. The postulates used to construct this conceptual model were justified using analysis of relevant physical and biogeochemical processes and feedbacks. In particular, the role of climate–ice sheet–carbon cycle feedback in shaping and globalization of glacial cycles is discussed. The reasons for the instability of the large northern ice sheets and the mechanisms of the Earth system escape from the “glacial trap” via a set of strongly nonlinear processes are presented. It is also shown that the transition from the 41 to the 100 kyr world about 1 million years ago can be explained by a gradual increase in the critical size of ice sheets, which in turn is related to the gradual removal of terrestrial sediments from the northern continents. The implications of this nonlinear paradigm for understanding Quaternary climate dynamics and the remaining knowledge gaps are finally discussed.</p

Sensitivity simulations with direct shortwave radiative forcing by aeolian dust during glacial cycles

Possible feedback effects between aeolian dust, climate and ice sheets are studied for the first time with an Earth system model of intermediate complexity over the late Pleistocene period. Correlations between climate and dust deposition records suggest that aeolian dust potentially plays an important role for the evolution of glacial cycles. Here climatic effects from the dust direct radiative forcing (DRF) caused by absorption and scattering of solar radiation are investigated. Key elements controlling the dust DRF are the atmospheric dust distribution and the absorption-scattering efficiency of dust aerosols. Effective physical parameters in the description of these elements are varied within uncertainty ranges known from available data and detailed model studies. Although the parameters can be reasonably constrained, the simulated dust DRF spans a~wide uncertainty range related to the strong nonlinearity of the Earth system. In our simulations, the dust DRF is highly localized. Medium-range parameters result in negative DRF of several watts per square metre in regions close to major dust sources and negligible values elsewhere. In the case of high absorption efficiency, the local dust DRF can reach positive values and the global mean DRF can be insignificantly small. In the case of low absorption efficiency, the dust DRF can produce a significant global cooling in glacial periods, which leads to a doubling of the maximum glacial ice volume relative to the case with small dust DRF. DRF-induced temperature and precipitation changes can either be attenuated or amplified through a feedback loop involving the dust cycle. The sensitivity experiments suggest that depending on dust optical parameters, dust DRF has the potential to either damp or reinforce glacial–interglacial climate changes

An efficient regional energy-moisture balance model for simulation of the Greenland Ice Sheet response to climate change

In order to explore the response of the Greenland ice sheet (GIS) to climate change on long (centennial to multi-millennial) time scales, a regional energy-moisture balance model has been developed. This model simulates seasonal variations of temperature and precipitation over Greenland and explicitly accounts for elevation and albedo feedbacks. From these fields, the annual mean surface temperature and surface mass balance can be determined and used to force an ice sheet model. The melt component of the surface mass balance is computed here using both a positive degree day approach and a more physically-based alternative that includes insolation and albedo explicitly. As a validation of the climate model, we first simulated temperature and precipitation over Greenland for the prescribed, present-day topography. Our simulated climatology compares well to observations and does not differ significantly from that of a simple parameterization used in many previous simulations. Furthermore, the calculated surface mass balance using both melt schemes falls within the range of recent regional climate model results. For a prescribed, ice-free state, the differences in simulated climatology between the regional energy-moisture balance model and the simple parameterization become significant, with our model showing much stronger summer warming. When coupled to a three-dimensional ice sheet model and initialized with present-day conditions, the two melt schemes both allow realistic simulations of the present-day GIS

Simulation of the last glacial cycle with a coupled climate ice-sheet model of intermediate complexity

A new version of the Earth system model of intermediate complexity, CLIMBER-2, which includes the three-dimensional polythermal ice-sheet model SICOPOLIS, is used to simulate the last glacial cycle forced by variations of the Earth's orbital parameters and atmospheric concentration of major greenhouse gases. The climate and ice-sheet components of the model are coupled bi-directionally through a physically-based surface energy and mass balance interface. The model accounts for the time-dependent effect of aeolian dust on planetary and snow albedo. The model successfully simulates the temporal and spatial dynamics of the major Northern Hemisphere (NH) ice sheets, including rapid glacial inception and strong asymmetry between the ice-sheet growth phase and glacial termination. Spatial extent and elevation of the ice sheets during the last glacial maximum agree reasonably well with palaeoclimate reconstructions. A suite of sensitivity experiments demonstrates that simulated ice-sheet evolution during the last glacial cycle is very sensitive to some parameters of the surface energy and mass-balance interface and dust module. The possibility of a considerable acceleration of the climate ice-sheet model is discussed

Simulation of the last glacial cycle with a coupled climate ice-sheet model of intermediate complexity

A new version of the Earth system model of intermediate complexity, CLIMBER-2, which includes the three-dimensional polythermal ice-sheet model SICOPOLIS, is used to simulate the last glacial cycle forced by variations of the Earth's orbital parameters and atmospheric concentration of major greenhouse gases. The climate and ice-sheet components of the model are coupled bi-directionally through a physically-based surface energy and mass balance interface. The model accounts for the time-dependent effect of aeolian dust on planetary and snow albedo. The model successfully simulates the temporal and spatial dynamics of the major Northern Hemisphere (NH) ice sheets, including rapid glacial inception and strong asymmetry between the ice-sheet growth phase and glacial termination. Spatial extent and elevation of the ice sheets during the last glacial maximum agree reasonably well with palaeoclimate reconstructions. A suite of sensitivity experiments demonstrates that simulated ice-sheet evolution during the last glacial cycle is very sensitive to some parameters of the surface energy and mass-balance interface and dust module. The possibility of a considerable acceleration of the climate ice-sheet model is discussed

Quantifying the effect of vegetation dynamics on the climate of the Last Glacial Maximum

International audienceThe importance of the biogeophysical atmosphere-vegetation feedback in comparison with the radiative effect of lower atmospheric CO2 concentrations and the presence of ice sheets at the last glacial maximum (LGM) is investigated with the climate system model CLIMBER-2. Equilibrium experiments reveal that most of the global cooling at the LGM (-5.1°C) relative to (natural) present-day conditions is caused by the introduction of ice sheets into the model (-3.0°C), followed by the effect of lower atmospheric CO2 levels at the LGM (-1.5°C), while a synergy between these two factors appears to be very small on global average. The biogeophysical effects of changes in vegetation cover are found to cool the global LGM climate by 0.6°C. The latter are most pronounced in the northern high latitudes, where the taiga-tundra feedback causes annually averaged temperature changes of up to -2.0°C, while the radiative effect of lower atmospheric CO2 in this region only produces a cooling of 1.5°C. Hence, in this region, the temperature changes caused by vegetation dynamics at the LGM exceed the cooling due to lower atmospheric CO2 concentrations